INTRINSIC CONSTRUCTION OF INVARIANT FUNCTIONS ON 

SIMPLE LIE ALGEBRAS 
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Abstract. An algorithm for constructing primitive adjoint-invariant func- 
tions on a complex simple Lie algebra is presented. The construction is in- 
trinsic in the sense that it does not resort to any representation. A primitive 
invariant function on the whole Lie algebra is obtained by lifting a coordinate 
function on a Kostant slice of the Lie algebra. Such an intrinsic construction of 
invariant functions is most useful for the bigger exceptional Lie algebras such 
as the E's. The Maple implementation of this algorithm is outlined at the end 
and will be applied to these exceptional Lie algebras in a future work. 



1. Introduction 

Let g be a complex simple Lie algebra of rank I with adjoint group G. We 
recall that G acts on g by the adjoint action, and therefore on the algebra 7(g) of 
polynomials on g by its contragredient, that is, 

(1.1) (g ■ P)(x) = P(Ad g -ix), geG, Pe 9(g), x e g. 
Let 

m = n&f 

be the algebra of polynomials on g invariant under the above action of G. A well- 
known theorem of Chevalley |Che55j asserts that 1(g) is polynomial algebra on I 
homogeneous polynomials 1% , ■ ■ ■ , // , that is, 

7( )=C[/ ll --- ,7,]. 

We will refer to the Jj's as primitive invariant functions on g. Write the degrees 

(1.2) degl j = d j , j = 1, 

We will assume that the Ij 's are ordered in the sense that 

d\ < c?2 < • • • < di . 

The numbers 

(1.3) rrij = dj - 1, j = 1, • • • ,1, 

are called the exponents of g. 

Although the choice of the Ij 's is not unique, the degrees dj and hence the ex- 
ponents m,j are intrinsic to g which constitute important invariants (see |Che52j ). 
Our main objective in this paper is to give an algorithm to explicitly and intrinsi- 
cally construct a set of primitive invariant functions. We note that our invariant 
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functions are defined on the whole Lie algebra g. We also comment that our con- 
struction is uniform, explains the pattern for the exponents, and does not resort to 
any representation. The author has implemented his algorithm on Maple. 

The traditional way of obtaining such invariant functions on g is extrinsic by 
employing a faithful representation. Usually the first fundamental representation of 
the Lie algebra is used because of its small dimension. The representation expresses 
a general Lie algebra element as a matrix, whose dimension is the same as that of 
the representation space, involving its coordinates. The sum of the principal minors 
of this matrix with dimension equal to a degree dj (|1.2[) of the Lie algebra is one 
such sought-after primitive invariant function. To this author, there are several 
drawbacks to this approach. First, this construction uses the a prior information 
of the degrees dj of g without being able to provide any deeper reason, and in the 
case of Dm = so(4n) where the degree 2n has multiplicity 2, a special formula 
is needed for the Pfaffian. Furthermore for bigger exceptional Lie algebras, their 
representations are hard to be made explicit, and the enormous cardinality of the 
principal minors of a big matrix prevents this method from being efficient. In 
particular, explicit forms of invariant functions on Eg are only known up to the 
second one of degree 8 [CP07] . 

In view of the above, an intrinsic and uniform method is clearly desirable. Intu- 
itively speaking, our algorithm uses the restriction of the adjoint representation on a 
principal 5(2 subalgebra in g, and we gain independence from other representations 
and furthermore computational efficiency in this way. Together with his collabora- 
tor, the author plans to apply his Maple implementation of this algorithm to other 
invariant functions on the E's, with the goal of at least getting some interesting 
restrictions of the invariant functions, that is, carrying out step (i) in Theorem ll.lll 
The restricted functions can still be very interesting and useful for the exceptional 
Lie algebras, since for example they are the functions most relevant to integrable 
systems such as Toda lattices [Kos79] and generic Toda flows |GS99j . 

The foundation for our construction is Kostant's profound studies Kos59,Kos63, 
IKos78] on invariant functions, which we now introduce. Fix a Cartan subalgebra 
i) C g. Let A be the corresponding root system, A± a choice of positive/negative 
roots, and ir — {ai, • • • , a{\ the positive simple roots. Let Q = t) ® QeA g a be 
the root space decomposition, with g a generated by a root vector e Q . For a £ A+, 
let H a = [e Q ,e_J, and we require a(H a ) — 2 for the choices of root vectors. For 
1 < i < I, the H ai form a basis of f). 

The height (or the order) o(a) of a root a £ A is defined as 

/ 1 
(1.4) o(a) = y^n», ifa = ^n l a i . 



This also induces a height gradation 

(1-5) g = 0gt, flfe = (J) Q and g = !)■ 

k o(a)=k 

For x 6 gk, we write o(x) — k by abusing the notation and call k the height of x. 
Let n = QeA+ g a = ©fe>o 3k be the maximal nilpotent subalgebra of q, b = f) ©n 
the Borel subalgebra, and N the unipotent subgroup of G corresponding to n. 
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Define 

l 

(1.6) e = ^-a i - 

i=l 

Let s be a complement of [e, q] in g, that is, 

(1.7) fl =s©[e, fl ]. 

Then by |Kos59, Kos63 , sen, and dim(s) = I is equal to the rank. We call s 
a Kostant slice, and let {sj}' =1 be a homogeneous basis of s with respect to the 
height gradation (|1.5p . 

The following theorem summarizes several important results of Kostant on in- 
variant functions. 

Theorem 1.8 (Kostant [Kos59, Kos63, Kos78j ). The heights of the Sj are corre- 
spondingly the exponents of the Lie algebra q. That is, if we order the {sj} so that 
o{si) < o(s2) < • • • < o(si), then 

(1.9) o(sj) = mj, l<j<l 

Furthermore, there is a sequence of isomorphisms through restrictions 

(1.10) I(Q) ^y(e + b) N A?(e + s), 

where 7(e + b) N is the algebra of polynomials on e+b invariant under the N action, 
and T(e + s) is the algebra of all polynomials on e + s. 

The following is our main result. 

Theorem 1.11. (i) There is an explicit algorithm for constructing the in- 

verse to r2 in (jl.lOj) : 

r^ 1 : 3>(e + s) -> 9(e + b) N . 

More precisely, let the £j be the coordinates of a general point 

i 

' ■ X • ' ' •'• 

j'=i 

Then there is an algorithm for constructing I primitive invariant function 
Ij defined on e + b of degree dj such that 
I 

(1.12) I j (e + J2&i)=Zj> l<i<*. 

i=l 

(ii) Furthermore, there is an explicit algorithm for constructing the inverse to 
T\ in (|1.10|) swc/i t/iat i/ie invariant functions are defined on the whole q. 

We present our basic setup and compution techniques in Section 2. In Section 
3, we present our algorithms in the proofs of the two parts of our main Theorem 
II. Ill together with several propositions. In part (i), our algorithm for constructing 
r^ 1 lifts the values of the invariant functions from the slice e + s to e + b and then 
to the Cartan subalgebra h. In Proposition I3.18| we prove the Weyl invariance 
of the resulted function on h. In part (ii), wc present the similar algorithm for 
constructing r^ 1 to define the invariant functions on g. The complexity of this 
step is much bigger than the previous step, and to a large extent accounts for the 
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difficulty in getting the invariant functions on the whole Lie algebra. In Section 4, 
we outline how the author has implemented the algorithm on the software Maple. 

Because of the special roles played by exceptional Lie algebras in mathematics 
and physics, we expect our algorithm and the explicit invariant functions it produces 
to have applications in a range of areas such as integrable systems and higher 
Casimir operators. 

Acknowledgment. The author thanks Ian Anderson for help with Maple and 
Luen-Chau Li for interest in this work. 



2. Setup and computation techniques 

Our later calculations rely on the following setup of Kostant |Kos63] in an essen- 
tial way. We have also extensively used such a setup and developed some techniques 
for computation in [LNllj . This paper is a further application of such techniques. 

To begin with, we recall that the polynomial algebra 7(g) can be identified with 
the symmetric algebra S = S(g*) on g* , the dual of g. On the other hand, we can 
associate to each i£ga differential operator d x on g, defined by 

(2-1) (&/)(») = 4 f{v + tx), /eC»(g). 

at t—o 

In this way we have a linear map i4fl r which can be extended to an isomorphism 
from the symmetric algebra = S(g) on g to the algebra of differential operators 
d with constant coefficients on g. From now onwards we will identify these two 
spaces. With this identification, we have a nondegenerate pairing between 5** and 
S given by 

(2.2) (d,f) = (df)(0), 9eS„ feS, 

where (9/)(0) denotes the value of the function df at G g. It is clear that both 
5* and S are graded from the tensor structure: 5* = (Bk>oS*, S — ®k>oS k , and 
Si pairs nontrivially only with S J . 

If / G S k and x G g, it follows from the Taylor expansion that 

( 2 -3) (W f ) =f{x) - 

It is clear that the adjoint action of G on g can be naturally extended to an 
action of G on 5*. On the other hand, S is a G- module as its contragredient by 
(jl.ip . (We denote the actions of G and later of g by a dot.) We have 

(2.4) (g-d, g-f) = (d,f), Vg G G, d G S„ f G S. 

By differentiation, S and S* become g- modules and the actions of g on both spaces 
are by derivations. Therefore we have the following properties: 

(2.5) [i, y] ■ d = x ■ (y ■ d) — y ■ (x ■ d), x,y G g, 3 6 5,, Lie alg horn 

(2.6) x ■ d y = d[ x y ], x, y G 0, adj action 

(2.7) x ■ (dS) = (x ■ d)5 + d(x -5), x G g, d, S G S*, derivatioin 

(2.8) x-d n = nd n ~ 1 (x-d), power rule 

(2.9) / G 1(g) x ■ f = 0, Vx G g. inv property 
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Since the pairing between 5* and S obeys (12. 4p . it follows from derivation that 
the g-actions satisfy 

(x-d,f) + (d,x- f) = 0, Vxeq, SeS», feS. 

This and f2~g]) imply that 

(2.10) (x-d,f) = 0, V/£l(g),i6j. 

There is a grading element so £ f) defined by the conditions that 

(2.11) oci(x ) = 1, V 1 < i < Z. 

By (|2.1ip . a(a;o) = o(a), and [xo,e Q ] = o(a)e a . Thus the graded subspaces from 
(|1.5|) are $k — {% & q\[xq, x] = kx}. This motivate the following definition of the 
weight structure of [Kos59j on S*. For each k £ Z, define 

(2.12) Sfc = {9 € S* | a; • d = kd}. 

Applying (]2 . 10[) to Xq g f) gives us the first vanishing result of Kostant [Kos63 . 
If d g 5 fc for fc ^ 0, then by (|2T2j) . d = ±x ■ d. Then in view of (I2TTU)) . we have 

(2.13) 5 g S k for fc ^ => (0, /) = for all / 6 7(g). 
Applying (12.10)) to a general x € i) gives us the following refined vanishing. 

Lemma 2.14 f |Kos78] . [LNTTl Lemma 3.1]). For all f g 7(g), pet), 



3; ncM>=° 



unless ^2 aeX TO " a = 0. 



In our later constructions, we will further exploit (12. 10)) by applying it to other 
elements in g. The most convenient formulation for us is the following "integration 
by parts" formula to "move things around." 

Lemma 2.15 f |LN11[ Lemma 6.2]). Let f g 7(g), then for all x, y g g, d G 5*, 

and m > 0, we have 

(2.i6) m*,^,/) = ^(sr +l (rs),/). 

For the reader's convenience, we repeat the brief proof. 
Proof. By using (|23]h (|2T8]) . flUTJ) and (j2~10| . we find that 

(dTd [x , y] dj) = -(d™(yd x )dj) 

1 -((y-9r +1 )5,/) 



771 + 1 

1 -(y ■ (d™ +1 d),f) + -^—(d x n+1 (y -d),f) 

1 777 +1 



777. 
1 



-{d™ + \yd),f). 



777 + 

□ 
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3. Constructive proof of the main theorem 

In this section, we present our algorithms in the forms of proofs to the two parts 
of our main Theorem II .111 and we present several supporting propositions. At the 
end of this section, we show how to assemble the information we obtain to get the 
concrete form of an invariant function. 

First we present a different basis of g, which is crucial to our inductive procedures 
later. Define 

(3.1) s*? = (ade)* s h 0<k<2m j , l<j<l. 

Then s) = [e^ 1 ] for k > 1. 

Lemma 3.2. { S J } i<^< f™ 3 * s a basis of g. 

Proof. This follows from (|1.7p and Kostant's work in Kos59]. □ 
The height of Sj is rrij — k by (|1.9|) and (|1.6p . and therefore 

(3.3) U := {s h j 'lol^l™' -1 is a basis of n = £f) g fc . 

fc>0 

We will denote a general element of It by u. Note that for such it's, either u G s or 
u e [e,g] in the decomposition (|1.7[) . (Actually the preimage of such a u g im(ad e ) 
is unique since ad e has no kernel in b |Kos59j .) 
Similarly we have that 

(3.4) W := {s^Tgif ^ * a basis of 9k . 

fc<-2 

We will denote a general element of W by w. Note that all such w's belong to [e, q] 

in sun . 

By the special nature of go = f) and g_i, this author will use their natural basis 
{i? ai }i<j</ and {e_ Q!i }i<j<; in the algorithms. 

Proof of Theorem l 1 . 1 1\ (i). Since we want to construct invariant functions, we will 
enforce the invariance property (|2.10j) or more explicitly Lemma l2.15l In this part, 
we will inductively show that this rule and the condition (|1.12|) determine the 
invariant function Ij on e + b. Then we prove that the constructed function when 
restricted on f) is invariant under the Weyl group in Propositions 13. 18l Proposition 
13.201 quickly shows that such functions are algebraically independent. 

In this proof, we will work with a fixed Ij, and we often write / for short. Also 
d = dj and m — rrij . 

Using (|2.3[) and the multinomial theorem, to describe / G T(e + b), we need to 
determine all the following polynomials 

(3-5) (dud b e d%, I) = (d u{l) d u{2) ■ ■ ■ d u{c) d b e d%, I) 

of degree a in p S h. Here U = (u(l), • ■ • ,u(cj) is a sequence with each u(i) 
(possibly repeating) from 11 in Q3.3p . The expression (I3.5P is nonzero only if 

a + b + c = d, 

(3 - 6) E°w i ))= & ' 

i=l 

by (|2.2p and (|2.13p . (This is related to the xo-grading by Kahzdan in |Kos78j .) 
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We run increasing induction on a and decreasing induction on b to define the 
polynomials in p.5[) . 

The first case is a = and b = d — 1 = m. Then we need to determine all 
the (9 U <9™, I) with o(u) — m. For u = Sj where j is our fixed index, applying the 
multinomial theorem and the vanishing results, we have 

L J (e + s 3 ) = j [ ({d S] +d t ) d ,I ) hy^M 

(3-7) =±(d Si d?J 3 ). by(EI3D 

Therefore the defining condition (|1.12l) in this case, Ij(e + Sj) = 1, implies 

(3.8) (d Sj d™,Ij)=m\. 

For u — Si G s with o(si) — m but i =/= j (hence the multiplicity of m as an 
exponent is at least 2 and this happens, among all the simple Lie algebras, only 
for I?2n and m = 2n — 1, n > 2), similarly to (|3.7p . the defining condition (|1 . 12|) in 
this case, Ij(e + Sj) = 0, implies 

(3.9) (d Si d™, I,) = 0, o( Sl ) = o( Sj ), i ± j. 

For u — [e,v] £ [e,fl], the vanishing property f|2 . 10[) . together with p.6[) and (|2.8I) . 

forces 

(& u d:\i) = (d [tM d:\i) = — 1 —{ V • di n+ \i) = o. 

m + 1 

Now we compute (|3.5p for a = and all b. (Such expressions for (|3.5[) are all 
numbers.) If in U, at least one w(i) £ [e,fl], say u(l) = [e, Ui], then by Lemma [2. 151 
we have, with U = (u(2), ■ ■ ■ ,m(c)), 

(dud^i) = (d u(1) d fJ d b e j) = (d [e , Vl] d d b e ,i) = ^l T (9 e h + 1 ( Vl . (da)), i 

(3-10) ! ,» 

= ^r7(z_/^ +1,9 «(2) • •• 9 [v 1 ,u(n)] ' ••d«(c),-0 
n=2 

where all the terms on the right have (6 + 1) e's, and hence are known from the 
induction hypothesis. We need to show compatibility when there are two u(i) e 
[e,g], and this is done in Proposition 13. 16l 
On the other hand, (|1.12j) implies 

(3.11) {dud b e , I)=0, if c > 2 and all the u(i) = s k G s. 



Now assume that the (|3.5j) have been computed when the degree in p is < a — 1 
with a > 1, and we compute it for degree a. 

We in effect use the fact that every p G f) is in [e, g] in (|1.7|) . Actually for 
p = J2\=iPi H ai, define 

(3.12) x p = '^2p l e ai , then - [e, x p ] = p. 

i=l 

Here be regarded a linear function in p with values in Si. 
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Then Lemma \2. 151 gives 

(dud b e d a p ,i) = -(dud b e d%-%, Xp] ,i) 



(3.13) 



1 

u — i22( d e +ld u(l) ■■•d[u{n),x p ) ■■■d u{c) d^- 1 ,I) 
n=l 

+ l T ^(d b e +1 d u d [p , Xp] d;- 2 1 i) 

^ I C 

^-jJ2J2Pi( d e + ld u(l)---d[u(n),e^]---du(c)dp^ 1 ,I) 



a-1 



1 n=l 
/ 



by (|3.12[) and hence 



which has degree 2 in p. The factors in the first sum have degrees a — 1 in p, and the 
factors in the second sum have degrees a — 2 in p. Using the induction hypothesis 
and p. 31) , all such factors can be expressed in the known cases of (|3.5[) . 
We can continue (|3.13[) all the way until we get a — d, where we have 

i 

(3.14) (d d p ,I) = (d-l)J2PiMp)(ded eoii d^- 2 ,I}. 

Then by (|2.3I) . our function I(p) on () is 

(3.15) J(p) = 

Proposotions 13.181 and [37201 below prove that the constructed this way 

are algebraically independent and invariant under the Weyl group when restricted 
to the Cartan subalgebra f). □ 

Proposition 3.16. There is compatibility when there are two choices for u G [e,g] 
in (j3~TU|) . 

Proof. Assume, for example, u{\) = [e, Vi],u(2) = [e, 1*2]. Then, with U' = 
(u(3), • • • ,u(c)), 

(dud b e) I) = (d u(1) d u{2) d w d b e) I) 

can be computed in two ways using v\ or w 2 in (|3. 10[) . The first answer A\ using 
Vi is, by (E7I1) and (l2~B . 

^i = ^t(^ +1 ( V i-(9„( 2 )%0)^) 

= ^ T ((5 e b+1 % 1 , tl(2)] ^,/) + (d b+1 d u{2) { Vl ■d u ,),I)) 
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and similarly for the second answer A 2 using v 2 . Therefore for the difference, we 
have 

(b + l)(A 1 -A 2 ) 

= (d b e +1 d {[viM2)Mv2Mm d w ,l) + (d^d^^ ■ du>),I) - {d b+1 d u{1) (v 2 ■ dur),I) 
= (dl^d^^du,, I) + (d^d^ivj ■ du,),l) - (d b e +1 d [etVl] (v 2 ■ du,),l) 

= ^ r2 ((d b e +2 {[v l M ■ du,),l) + (d b+2 (v 2 ■ ( Vl ■ du,)),l) 

-(^+ 2 ( vl .(« 2 .^)),l)) 

= 0, 

where the first term in the second equality uses the Jacobi identity 

[vi,u{2)} - [v 2 ,u(l)} = [v!, [e,v 2 ]} + [[e,vi],v 2 ] = [e, [vi,v 2 ]}, 

the third equality uses Lemma |2 . 1 51 again, and the last identity uses the Lie algebra 
homomorphism property (|2.5|) (with its root in the Jacobi identity). □ 

Remark 3.17. Furthermore when a > 1, if in U there exists a u G [e, fl], there is 
an alternative approach similar to (|3.10p . which is compatible with (|3.13j) . by the 
same reason as above. 

Proposition 3.18. The I(p) on f) defined in (|3.15[) is invariant under the Weyl 
group W . 

Proof. Since W is generated on t) by simple reflections through the hyperplanes 
defined by en = for 1 < i < I , we only need to prove the invariance of the function 
I(p) under 7^ for any i. Fix an i and we omit it from the notation. 

We use an orthogonal basis of f) with the first vector being H a = H Ui . Then 
we write 

(3.19) p = xH a +Y, Y±H a ^a(Y) = 0. 

In this orthogonal basis, the reflection is just the transformation x H> — x, and 
we only need to prove the (dp, I) in p,14[) is a function of x 2 . For that purpose we 
run decreasing induction on k and I to prove that in general the 

D(k,l) := {d k ea d k e _ a d l Y d d p - 2k -\l) 

are functions of x 2 , with (3^,1) = D(0, 0). 

When2fc+Z > d, D(k,l) = by $TZS). When2/c+Z = d, D(k,l) = {d k a d k _ a d l Y , I) 
are constants with respect to x, since Y doesn't involve x. 
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Now by Lemma \2. 151 and using (|3.19|) . we have, for d — 2k — I > 1, 
D(k,l) = (dldtAd^-'-^xd^ +dy),I) 

= x(d Ha d k e jt a d l Y d*- 2k - l -\l) + (dldtj^d*- 2 '- 1 - 1 ,!) 
= x(d [e ^ e _ a] d k ea d k e _ a d l Yd d p ~ 2k ^\l) + D(kJ + 1) 

1 l( d e-l( e « ■ (dt^- 2k - l - 1 )),l)+D(k,l + l) 

xa(p)(5 e fc + 1 5 e fc + 1 ^^- 2fc - i - 2 , /> + D{k, I + 1) 



k 

d-2k-l-l 



k + 1 

■ / - 1 ^ „ 

x 2 D(k + l,l) + D(k,l+1), 



2{d-2k-l- 1) 2 



since 



-ea.^*,^- 3 *-'- 1 ) 

= ^^-^^ej^"^-'- 1 + (d- 2k - I - l)d k a d l Y df- 2k - l - 2 d { ^ ea] 

= (d-2k-l- l)a{p)d k + 1 d Y d^- 2k - 1 - 2 

due to that [Y,e a ] = a(Y)e a = 0, [p, e a ] = a(p)e a , and a(p) — 2x by (|3.19[) . 
Therefore the appearance of x in D(k, I) is always through an x 2 entry. □ 

Proposition 3.20. The are algebraically independent. 

Proof. This is clear from our defining condition (11.121) . since the Ij restrict to the 
coordinates £j on the slice e + s. □ 

Remark 3.21. In a sense, the above algorithm is the reversion of the procedures in 
[LN111 §6]. Here we start with a high root vector and push the function down to 
f). In [LNlTj we derived information about higher and higher root vectors starting 
from some knowledge on t). The direction here is more delicate. 

Proof of Theorem A 1 . 1 1\ (ii). The further lifting of the invariant function / = Ij to 
the whole Lie algebra q involves considering all such terms 

(3.22) (d w d b t d;du,I) 

where dw = 9 w m ■ ■ ■ d w ^ with each w(i) from W in (|3.4|) . and du is the same as 
in (|3.5[) . The absolute value of the total weight of dw is 



(3.23) -o(W) = -^2o( Wi ). 

i=l 

Similarly to (13.61) . (|3.22p is nonzero only if 

(3 + a + b + c = d, 

(3 ' 24) -o(W) + b = J2o(u(i)). 

i=l 

Note that we do not need any vectors with heights or —1 in p.22j) . since all 
such vectors are accounted for by the <9° and d b terms using the following Lemma. 
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Lemma 3.25. If a > and one copy of the p € f) is replaced by H ai , then 

(3.26) (dwd^d^dn^ d u,I) = \ d Hai {(dwd^du, I)) . 

Combining the U and W in (|3.22|) together and changing them back to the stan- 
dard basis using root vectors, let V — {e aii , • • • , e ai } be a sequence of (possibly 

repeating) root vectors from fe>o 0fc©0 fc <_ 2 0k- Let ay = J2]=i a ij = Si=i n i a i 
be the sum of the roots and express it in terms of the simple roots. Assume that 
ni > for 1 < i < I. Let dy = d Ca . ■■•d ea . , and let o{V) ~ X)!=i n « be the 
grading of dy in as defined in (|2.12[) . Let I £ /(g) be an invariant function 

with degree d > 7 + o{V). We have 

(3.27) (d°^d v d d -^ v \l) = n °W l n / ( J] de"!, ) dyd^-°W , A , 

* i—l ' 

where pel). 

// o(V ) > and one copy of e is replaced by e_ Qi , then 

(3.28) (^^'"^e^^y^ 7 " ^',/) = -^{d^dyd^-^J). 

Proof of Lemma \3. 251 Let p = Yl%=i ViH ai G f). The on the right hand side of 
(|3.26l) stands for . All these formulas are easy consequences of the multinomial 
theorem for the powers of d p and d e , and in (|3.27[) and (|3.28[) we also need to use 
Lemma \2. 141 We omit the details. □ 

Now returning to the main proof. We run induction on the lexicographical order 
of the pair (/?, —o(W)), with —o(W) defined in (|3.23[) . to determine such terms in 
(|3.22l) . Since w(l) <E [e,g] from (|3~4l) . assume 

(3.29) w(l) = [e,v 1 }. 

Then v x 6 © fc <_ 1 fl fe . With W = (w(2), ■ ■ ■ ,w(fi)), Lemma EE] gives 
{d w d b e d a p du,I) = (d^d^d^duj) 



1 P 



/3 

(2) • • •<%>!, «>(m)] • • • d w (p)d b t +1 dpdu , I) 



m=2 

3-30 a , , 

1 c 

n=l 

Here all the summands in the first sum can be expressed by (|3.22[) with /3 — 1 
elements from W in p.4[) . The summands in the second and the third sums can 
either be expressed by (|3.22p with f3 — 1 elements from W if the heights of [vi , p] 
or [vi,u(n)] are > —1 by Lemma 13.251 or with j3 elements from W otherwise. But 
the new — o([«i,p]) or — o([vi, u(n)j) is strictly less than the old — o(w(l)), since the 
height of vi is one bigger than that of w(l) in view of (|3.29j) . and p and u(n) have 
nonnegative heights. Therefore the new total — o(W) is strictly less than the old 
one. 
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Therefore through this hierarchy of induction hypothesis, all the terms on the 
right are known. 

We note that the outcome of (|3.30|) does not depend on the choice of v\ in 
(|3.29p . which may not be unique. Say v[ = Vx + vo with [e,vo] — 0. Then the 
outcome of (I3.30P is linear in v\ and vq, and the terms for combine to give 
{dit^dyyd^dpdu, I) = by tracing the identity backward. □ 

After all these coefficients in (|3.5[) and (I3.22p . as functions on p £ f), are calcu- 
lated, we can assembel our function as follows. Let 

x = 2J z i w i + e + p + 2J Vj u j 

be an element in g with the Wi from W in (|3.4p . the Uj from XL in (|3.3|) . pel), and 
the Zi and yj as coefficients. Then we get I(x) by (|2.3[) . the multinomial theorem, 
and the coefficients (|3.5[1 and (|3.22p . 

If we change the basis back to the usual root vectors, then we get I(x) for 
x = p + e + Yl (a)^-i x a e a- Using (|3.27l) . we can further spell out the dependence 
on the e-OLi ■ At the end, we obtain the function I{x) expressed in the coordinates 
of a general element in q: 

i 

X — ^ Pi H a i -\- ^ X a G a . 
i=l QgA 

See the end of Section 4 for an example. 

4. Implementation of the algorithm on Maple 

It turns out that our algorithm is very ready for implementation on Maple, espe- 
cially using the LieAlgebras package under Maple written by Prof. Ian Anderson. 
One particularly useful feature is that we can do the change of basis in Lemma 
13.21 easily. In this section we outline the steps that we run on Maple. This author 
has written a Maple program containing all the implementations. He is currently 
working with his collaborator to apply the program to the exceptional Lie algebras 
of the E type. 

First we setup the Lie algebra g in Maple, preferably under a Chevalley basis. 
This is achieved by either giving a list of root vectors as matrices if one knows a 
representation, or by inputting the non-zero structure equations. For the excep- 
tional Lie algebras of the E type, the structure constants for a Chevalley basis are 
available, for example following t Vav04 a . After the DGsetup command in Maple, we 
have the basis vectors labeled by ei for 1 < i < dimg and the Lie bracket operation 
available. The corresponding roots for the ei are readily calculated in terms of 
the simple roots, and so are the heights (jl.4p . This author has made the choice of 
ordering the root vectors according to the heights from to the maximal height, 
and then from —1 to the negative maximal height. 

Then we make the change of basis from the {e±} to the basis in Lemma [3.21 
Here we need to input the homogeneous slice basis {sj} (|1.7p . In practice, it 
is sufficient to choose root vectors at each height which appears as an exponent 
according to its multiplicity. The e in (| 1 .6(1 is readily calculated, and by inductively 
calculating (|3.ip . we get a new basis. We order these basis vectors according to 
their homogenous heights as before, and set up the Lie algebra again in this new 
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basis. This is particularly easy using the LieAlgebras package, and in particular 
we have the transition matrix M immediately ready. It is this matrix M that allows 
us to go back and forth between the two bases easily. We also record where we have 
the relation u = [e,v] in (|3.1I) . and construct a procedure fv to find a preimage of 
an element in li (|3.3[) or W (|3 .4[) under ad e . 

The first big step for the program is to find all the terms (|3.5|) and more generally 
()3.22|) which are possibly nonzero. In this implementation, we record each term 
(13.221) as a polynomial using the W and U information, and we further record its b 
and a values, which are dictated by p.24[) . 

For this we use the vanishing condition (|2.13p . since our new basis vectors from 
II in (|3.3p and W in (|3.4|) are still homogeneous with respect to the height grad- 
ing (|1.5[) . First we compute the possible height vectors for W and U, and then 
we readily get all the monomials corresponding to W and U with these height 
structures. 

Let mi (jl.3l) be the maximal height of the roots. Using the absolute values of 
the heights for dw for simplicity, we are looking for partitions \oY , 02 ■, • ■ ■ , of 
the total height —o(W) in (I3.23|) into f3 parts. Here we require that 2 < of 7 < mi 
for 1 < i < (3, For the heights of djj, we are looking for partitions [o^ , • ■ ■ , o^] of 
the total height — o(W) + b determined by (|3 . 24[) into c parts with 1 < d( < mi. 

All such termss can be generated using modified partition functions by running 
through some loops. See Appendix A. We store the results in several lists. The 
list ttms stands for top-terms and corresponds to /3 = and a = in (|3.22|) . The 
list ptms stands for p-terms with p £ I) and corresponds to /J = and a > in 
(|3.22p . The list ntms stands for negative-terms and corresponds to (3 > in (|3.22|) . 
In Appendix A, we have paid attention so that the orders of the lists follow our 
inductive procedures in our proofs in Section 3. 

The partition function grows very fast, and this corresponds to the fact that 
the number ~^7^Erj! °f terms in a degree d polynomial in n variables grows very 
fast with d and n. Therefore it takes a huge computational time to complete this 
generating procedure in Appendix A when the degree d and the dimension of the 
Lie algebra n get big. If we are only interested in getting a function in CP(e + b) as 
in part (i) of Theorem 1 then we can let /3 — 0. The corresponding procedure 
is much shorter, and we still get interesting functions in y(e + b) N . 

The terms in (|3 . 8|) , (|3.9I) and (13. lip are characterized by the fact all the factors in 
djj do not come from [e, g]. We easily find them in our list ttms using the procedure 
fv. Although we have specified their values in (|3.8p . (I3.9P and (|3.1ip . we can 
actually define the corresponding values as arbitrary constants, which corresponds 
to getting an algebraic combination of the primitive invariant functions. We store 
the values in a list valuedata. 

Then we run through our inductive algorithms for the lists ttms, ptms and ntms 
using the formulas (|3~T0l) . (|3~13)) and (EHO)) . Note that for the ((330)) terms, we 
need to take into considerations Lemma 13.251 Each time we compute a term using 
the previous valuedata information, and we add the result to the list valuedata. 
These formulas are quite ready to be coded up, as this author has done. There is 
an extra term which consists of purely Cartan terms, and it is calculated by (13.141) . 

After we get all the terms stored in valuedata, we can put the function together, 
change back the basis, and spell out the dependence on the e_ Qi as explained at 
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the end of Section 3. Our end result is an invariant function I(x) for x e in its 
standard form. 

Now we illustrate our program by the degree 6 invariant function on 532. 
We use the basis of 32 as made explicit in the Appendix of |BFO + 90] . We setup 
our 02 with basis 



ei 


= H ai 


e2 = 


Ha 2 




e3 


— e ai 


e 4 = 




es e Ql -|_ a2 


e-% 


— 62ai+CK2 


e7 = 




e8 = e3 Ql +2a2 


eg 


— 6_ ai 


eio = 


- e_ a2 





^12 6_2ai— Q2 ^13 6_3 Q1 _ a2 e^4 e_3 Q1 _2a 2 

We choose our slice elements in (jl.7l) to be si = 64, s 2 = eg. Let e = eg + eio, 
and we do the change of basis in Lemma l3.2l Order the new basis according to our 
above convention, and denote them by ■ 

The degree d is set to be 6. A nonzero term from (|3.22[) . for example, 

(4.1) (dj 11 df s d e dp,I) is recorded by yfiUs with b — 1 and a = 2. 

We generate the possible nonzero terms by the procedure in Appendix A. There 
are 8 terms in ttms, 10 in ptms with one extra purely on the Cartan subalgebra 
(|3.14[) , and 535 in ntms. 

Among the terms in ttms, the inputs are calculated to be ys and corresponding 
to (df s d^,I) and (<9§t?f,7). We require that their values be 5! and by (I3.8|) and 
p. lip , that is, we start our valuedata = [[y 8 ,240], 0]]. Running through our 
program we get all the values for the lists ttms, ptms and ntms. 

For example, our algorithm and program compute that 

(d 2 fii d f8 d e d 2 p ,I) = 3136 (3p 2 -Pa)(3pi - 2ft,). 

Then by ()2.3j) and the multinomial theorem, I(p + e + ^2 i — 3 Uifi + ^j— 11 Vjfj) has 
one term 

^(dj^df^dl^ys = 784y2 l2/8 ( 3p2 - p 2 )(3 Pl - 2p 2 ). 

Now change the basis back to the root vectors ej, I(p + e + ^*_ 3 XiBi + X}jiii x j e j) 
has one term 

£11 a^8 (3^2 - P2)(3pi - 2p 2 ), 
since the relevant change of basis is /§ = es and fn = 28 en. Finally spelling out 
the dependence of the function on the e_ Qi by (|3.27|) . we see that IQ2i—iPiH ai + 
X)i=i + Sjin ^i) has one term 

x^xgx^Zpi -P2)(3pi - 2p 2 ), 

since the total root of d 2 lx d e& is 2{—oi\ — 012) + (3«i + 2a?2) = ai, and e_ Ql = eg in 
our ordering. 

A computer program can do all the above to all the terms in our valuedata, 
and eventually we obtain our function I(x) for a general 1 G 92- It turns out to 
be one quarter of the sum of principal minors of dimension 6 of the corresponding 
matrix representation of g2- The whole procedure takes about 10 seconds on a 
usual laptop. 
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With his collaborator, the author will apply this program to the E's, at least 
aiming to get the invariant functions on e + b. He believes that this program has 
computational advantage, since the corresponding functions are doable on a usual 
laptop. For example, he can compute the degree 12 invariant function on e + fa of 
Eq, while the usual method of taking the sum of principal minors can not be done 
on a usual laptop even in the degree 5 case when using the adjoint representation. 

Appendix A. Generate and order terms 

We explain more the Maple codes to generate all the possibly nonzero terms in 
(|3.22[) . Here d is the degree of the invariant function, and rax = mi (|1.3|) is the 
maximal height of roots. The author has actually worked out in more detail the 
ranges for all the indices. Here floor and ceil are the standard integer functions, 
ow = — o(W) is the absolute value of the total height of W (|3.23j) . The functions 
negpart (x,y ,z) and pospart(x,y,z) are the author's modified partition func- 
tions, implemented elsewhere, to partition x in to z parts with each part ranging 
from 2 for negpart (or 1 for pospart) to y. nwts and pwts are the absolute values 
of the negative and positive height vectors for dw and djj in (|3 . 22[) . Also allmul 
is the author's procedure, implemented elsewhere, to come up with the terms such 
as (|4.1|) in (|3.22[) from the height vectors. 

ttms: = [ ] : ptms: = [ ] : ntms: = [ ] : 

for beta from to floor (d*mx/ (mx+2) ) do 

for ow from 2*beta to min(beta, d-beta)*mx do 

nwts := negpart(ow, mx, beta): 

for a from to d-beta-1 do 

for b from max(0, ceil ( ( (d-beta-a) -ow) /2) ) to 
floor ( ( (d-beta-a) *mx-ow) / (mx+1) ) do 
pwts := pospart (ow+b, mx, d-beta-a-b) : 
if beta=0 then 
if a=0 then 

ttms := [seq(seq([a, b, f ] , f in allmul(ls)), Is in pwts), op(ttms)] : 
else 

ptms := [op(ptms), seq(seq([a, b, f ] , f in allmul(ls)), Is in pwts)]: 

fi: 

else 

ntms := [op(ntms), seq(seq(seq(seq( [beta, nt , f*g, a, b] , 

g in allmul(ls)), f in allmul (-nt) ) , Is in pwts), nt in nwts)]: 

fi: 

od: od: od: od: 
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